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ABSTRACT 


The  original  coherent  lightwave  systems  were  expected  to  offer 
significant  performance  gains  relative  to  standard  direct  detection 
systems.  This  expectation  has  not  been  realized  due  to  the  effects  of 
laser  phase  noise.  The  laser  phase  noise  process  results  in  the 
integration  of  a  random  variable  that  transitions  over  the  integration 
period  from  a  Gaussian  distribution  to  a  uniform  distribution.  The  use 
of  convolutional  coding  effectively  replaces  a  single  bit  time,  with  its 
mostly  noncoherently  integrating  latter  portion,  by  several  more 
coherently  integrating  bits.  This  primary  bit-time  effect  comes  in 
addition  to  the  normal  coding  effect  of  efficiently  trading  bandwidth  for 
error  performance.  The  improvement  in  performance  brought  about 
from  coding  may  enable  coherent  systems  to  live  up  to  previous 
expectations. 

The  contributions  of  this  thesis  include  the  visualization  of  the 
phase  noise  process,  the  efficient  computation  of  the  laser  phase  noise 
power  factor  probability  density  function,  and  the  computation  of 
performance  curves  for  uncoded  and  coded  systems.  Additional  sections 
on  coherent  lightwave  systems  and  coding  provide  tutorial  information. 
A  potential  military  application  is  discussed,  along  with  practical 
implementation  issues. 
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I.  INTRODUCTION 


Optical  communications  links  have  become  the  preferred  channels 
for  high  data  rate  applications  for  many  reasons.  A  single  fiber  can 
support  the  bandwidth  of  hundreds  of  coaxial  cables;  fibers  are  immune 
from  radio  frequency  interference;  damage  to  fiber  cables  is  easily 
located;  fiber  is  lower  in  weight  than  metal  cables.  As  fiber  cables  have 
become  more  rugged,  their  use  in  military  systems  has  become  more 
practical.  Fiber  is  now  being  proposed  for  guided  weapon  data  links, 
battlefield  communication  networks,  and  sensor  systems. 

The  current  commercial  state  of  the  art  is  direct  detection.  Direct 
detection  lightwave  systems  are  similar  in  concept  to  A.M.  radio,  in  that 
the  process  of  detecting  the  signal  makes  no  use  of  the  signal  phase.  In 
essence,  the  signals  consist  of  pulses  of  noise  energy,  with  a  bit  "1" 
represented  by  a  pulse  of  energy,  and  a  bit  "0"  represented  by  no  signal. 
The  preferred  channel  of  propagation  is  single  mode  fiber.  Fibers  are 
cylindrical  waveguides  subject  to  the  same  physical  laws  as  microwave 
waveguides.  The  term  "single-mode"  refers  to  the  support  of  only  one 
propagating  mode.  This  eliminates  dispersion  due  to  intermodal  mixing. 
Fibers  are  subject  to  attenuation,  which  is  logarithmic  with  distance. 
Fibers  also  suffer  from  chromatic  dispersion,  where  the  spectral 
components  of  the  source's  non  zero  linewidth  travel  at  different 
velocities  through  the  waveguide  material.  Chromatic  dispersion  effects 
are  reduced  through  the  use  of  laser  diodes  with  small  linewidths. 
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Attenuation  is  reduced  through  the  formulation  of  fiber  materials  and 
strict  elimination  of  water  and  OH'  molecules  from  the  fiber  during 
manufacturing. 

Improved  performance  is  theoretically  possible  using  coherent 
systems.  In  a  coherent  system  receiver,  the  source  light  signal  is 
optically  mixed  with  a  local  oscillator  optical  signal,  resulting  in  an 
electronic  intermediate  frequency  signal.  The  electronic  EF  signal  is 
processed  the  same  way  as  in  a  microwave  coherent  communication 
system.  The  source  and  local  oscillator  lasers  must  be  carefully 
controlled  to  maintain  an  EF  signal  compatible  with  the  demodulator. 
This  makes  coherent  systems  more  complicated  than  direct  detection 
systems.  Coherent  receivers  show  significant  sensitivity  improvements 
(up  to  20  dB)  as  compared  to  standard  direct  detection  systems. 

A.  THE  PROBLEM,  HISTORY,  AND  THIS  THESIS 

The  significant  barrier  to  the  use  of  coherent  systems  is  the 
presence  of  phase  noise  in  the  source  laser  and  in  the  receiver  local 
oscillator  laser.  The  phase  noise  results  in  source  linewidfh  broadening. 
The  nature  of  phase  noise  is  such  that  its  effects  decrease  with 
increasing  data  rates.  The  commercial  implementation  of  coherent 
systems  awaits  the  development  of  laser  sources  with  small  linewidth 
and  rapid  modulation  capability. 

The  effects  of  phase  noise  became  known  after  early  experiments  in 
coherent  systems  resulted  in  dismal  performance.  One  of  the  first 
comprehensive  papers  on  the  nature  of  phase  noise  was  published  by 
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Salz  in  1985  (Salz,  1985).  The  significant  obstacle  to  quantifying  the 
effects  of  phase  noise  is  the  determination  of  the  probability  density 
function  for  the  phase  noise  process.  Approximations  to  the  density 
were  put  forth  in  a  series  of  papers  by  Foschini,  Greenstein  and 
Vannucci,  (Foschini,  Greenstein,  and  Vannucci,  1988;  and  Foschini  and 
Vannucci,  1988).  Azizoglu  and  Flumblet  (Azizoglu  and  Humblet,  1991) 
refined  the  earlier  treatment,  and  demonstrated  how  multisampling 
receivers  can  combat  the  effects  of  phase  noise.  Since  multisampling  is 
equivalent  to  diversity,  and  since  diversity  is  equivalent  to  simple 
repetition  coding,  it  makes  sense  that  the  application  of  more  efficient 
forms  of  coding  can  provide  greater  benefits  than  diversify. 

This  thesis  is  about  the  application  of  convolutional  coding  to 
binary  frequency  shift  keying  (2FSK)  coherent  lightwave  systems  with 
noncoherent  detection.  It  is  shown  that  coding  offers  significant 
improvements  in  error  performance  for  a  given  data  rate. 

B.  CHOICES,  TRADEOFFS,  AND  COSTS 

The  choice  of  noncoherent  detection  is  driven  by  the  more  relaxed 
linewidth  requirements  for  this  type  of  system.  With  the  choice  of 
noncoherent  detection,  the  possible  modulation  formats  are  on-off 
keying  (OOK),  and  frequency  shift  keying  (FSK).  The  error  performance 
of  these  two  systems  referenced  to  average  bit  energy  to  noise  ratios  are 
the  same,  but  the  peak  power  required  in  the  FSK  system  to  achieve 
equal  average  power  is  3  dB  less  than  that  for  OOK.  The  type  of  coding 
considered  is  hard  decision  convolutional  codes.  While  soft  decision 


decoding  offers  several  dB  improvement,  it  is  more  complicated  and 
requires  a  decoder  that  is  more  difficult  to  implement  at  high  signal 
rates. 

The  pha^nioise  effect  decreases  with  data  rate,  and  it  will  be  shown 
that  for  uncodeu  systems,  the  drop  in  bit  energy  from  an  increase  in 
data  rate  is  generally  more  than  offset  by  the  drop  in  phase  noise  effect. 
Practical  systems  are  limited  in  signal  rate;  it  is  not  always  possible  to 
increase  coding  gain  without  encountering  this  signal  rate  limit.  This 
means  that  there  will  be  tradeoffs  between  increasing  data  rate  and 
employing  coding.  The  requirements  placed  on  the  system  must  be 
stated  in  terms  of  required  data  rate,  limits  on  signal  rate,  and  error 
performance.  In  many  circumstances,  employing  coding  will  be 
beneficial,  in  others,  it  may  not. 

The  tradeoff  options  mount  when  the  ability  to  frequency  multiplex 
coherent  systems  within  common  channels  is  considered.  Now,  several 
lower  data  rate  coded  channels  may  be  multiplexed,  providing  the  data 
rate  of  a  single  uncoded  channel,  at  lower  total  power  (but  greater 
complexity).  The  error  performance  curves  presented  in  this  thesis  can 
be  utilized  to  conduct  tradeoff  studies  between  single  channel  and 
multiplexed  systems. 

The  cost  of  implementing  coding  may  be  very  small.  For 
convolutional  coding,  the  encoders  are  very  simple.  The  decoders  are 
somewhat  more  complicated.  The  Viterbi  decoders  are  well  suited  to 
implementation  as  application  specific  integrated  circuits  (ASICs). 

Several  companies  advertise  decoders  for  low  data  rates  (10  Mbps)  for 
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$25  per  chip.  The  fastest  available  Viterbi  decoder  is  a  hard  decision 
system  capable  of  operating  at  signal  rates  of  200  Mbps  (NASA,  1993). 

The  capabilities  of  decoders  will  increase  as  logic  speeds  increase  and 
parallel  architectures  are  exploited. 

C.  THESIS  OUTLINE 

The  outline  of  the  thesis  is  as  follows.  The  first  s:  :tion  provides 
background  information  of  a  tutorial  nature  on  optical  communication 
systems,  and  the  phase  noise  process.  The  next  section  develops  the 
effects  of  phase  noise  on  uncoded  2FSK  systems  with  noncoherent 
detection.  New  contributions  in  this  area  include  the  visualization  of  the 
phase  noise  process,  and  the  efficient  computation  of  the  probability 
density  function  from  its  approximate  moments.  The  next  section 
begins  with  tutorial  material  on  convolutional  coding,  and  then  develops 
the  error  performance  for  coded  systems  using  several  simple  codes. 

The  concluding  section  considers  a  hypothetical  coded  coherent  system 
and  compares  it  to  a  conventional  direct  detection  implementation.  An 
appendix  contains  computer  codes  of  the  probability  density  function 
(pdf)  computation  and  code  transfer  function  determination  for  the 
Mathematica  computer  algebra  system. 
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II.  OPTIC  \L  COMMUNICATION  SYSTEMS 


Lightwave  systems  may  be  broadly  classified  as  noncoherent 
system^  or  coherent  systems. 

A.  NONCOHERENT  SYSTEMS 

In  a  noncoherent  system  (Figure  2.1),  a  source  laser  is  intensity 
modulated,  the  resulting  optical  signal  is  sent  over  multimode  or  single 
inode  fiber  cables,  or  propagated  through  free  space.  Fiber  optic  cables 
are  cylindrical  waveguides,  and  thus  light  propagates  under  the  same 
physical  laws  as  apply  to  microwave  waveguides,  i.e.,  the  waveguide  can 
propagate  either  single  or  multiple  modes.  Multiple  mode  fibers  are 
affected  by  the  interactions  between  the  modes,  causing  modal 
dispersion,  which  widens  the  spectrum  of  transmitted  light.  The  longer 
the  cable,  the  more  severe  the  modal  dispersion  becomes,  leading  to 
performance  degradation.  Obviously,  a  single  mode  fiber  does  not  suffer 
from  modal  dispersion.  Another  form  of  dispersion  that  all  fibers  suffer 
from  is  chromatic  dispersion,  where  different  frequencies  travel  at 
different  velocities  due  to  a  non  constant  index  of  refraction.  The 
solution  to  this  problem  is  to  use  a  source  with  as  small  a  linewidth  as 
possible. 
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Figure  2.1  Noncoherent  Lightwave  System 


When  employing  wide  linewidth  sources  (like  LED's)  over  multimode 
fibers,  the  maximum  link  distance  is  on  the  order  of  10  Km  (Keiser, 

1991).  At  the  receiver,  the  intensity  modulated  signal  power  is  directly 
detected  by  a  photodiode,  generating  a  voltage  proportional  to  received 
optical  power  plus  noise.  The  similarity  between  noncoherent  lightwave 
systems  and  a  crystal  A.M.  radio  are  now  clear.  Advantages  of 
noncoherent  systems  include  simplicity,  and  the  ability  to  use  relatively 
noncoherent  light  sources  such  as  light  emitting  diodes,  as  opposed  to 
more  coherent,  more  complicated,  and  more  costly  laser  diodes.  The 
disadvantages  of  noncoherent  systems  are  difficulties  in  multiplexing 
multiple  signals  on  multiple  lightwave  frequencies,  and  inferior  error 
and  inferior  repeaterless  range  compared  to  coherent  systems. 
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B.  COHERENT  SYSTEMS 


Coherent  lightwave  systems  (Figure  2.2)  employ  a  spectrally  very 
pure  source  laser,  and  single  mode  propagation  channels.  The  spectral 
widening  through  a  multimode  fiber  would  be  intolerable  in  a  coherent 
system,  so  current  systems  employ  single  mode  fibers.  Free  space,  low 
dispersion  propagation  is  also  possible.  At  the  receiver,  the  incoming 
laser  light  is  first  mixed  with  a  local  oscillator  laser  (Figure  2.3).  When 
the  mixing  is  to  baseband,  the  process  is  termed  homodyne  detection. 
When  the  mixing  is  to  an  intermediate  frequency,  it  is  termed 
heterodyne  detection.  The  detection  and  decision  process  is  now  carried 
out  electronically,  exactly  as  in  a  coherent  microwave  communication 
system.  Coherent  systems  show  up  to  20  dB  better  receiver  sensitivity 
compared  to  noncoherent  systems  (Keiser,  1991),  and  also  a  high  degree 
of  frequency  selectivity.  The  repeaterless  range  of  coherent  systems  is 
on  the  order  of  100  Km.  Figure  2.4  shows  how  a  large  number  of 
channels  (assuming  10  GHz  channel  spacing)  can  theoretically  be  fit  into 
the  windows  of  low  fiber  attenuation  available  with  current  optical 
fibers.  Further  information  on  fiber  optic  systems  can  be  found  in 
(Keiser,  1991). 

The  nature  of  coherent  systems  demands  close  control  over  the 
frequencies  of  source  and  local  lasers,  which  adds  complexity  to 
coherent  systems.  The  frequency  of  current  laser  diodes  is  affected  by 
the  laser  temperature,  and  the  drive  current.  Modulating  the  drive 
current  is  one  method  to  impart  frequency  modulation  on  the  source 
laser.  Typical  sensitivities  of  a  laser  are  10-20  GHz/°C  temperature  and 
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1-5  GHz/mA  drive  current.  Since  the  intermediate  frequency  is  usually 
m  the  MHz  range,  slight  deviations  in  temperature  or  drive  current  in 
source  or  local  lasers  can  cause  trouble. 


Figure  2.2  Coherent  Heterodyne  Lightwave  System 


Figure  2.3  Heterodyne  Optical  Receiver 


C.  OPTICAL  RECEIVERS 

A  heterodyne  optical  receiver  is  shown  in  Figure  2.3.  In  this 
receiver,  the  incoming  source  signal  is  projected  onto  the  photodetector 


f) 


(surface.  A  local  oscillator  laser  is  optically  mixed  via  a  mirror  onto  the 
photodetector.  Denote  the  source  signal  field  by 

Er  =Ar(t)cos{2xfrt  +  er(t)]  (2.1) 

where  fr  is  the  carrier  frequency  of  the  source  laser,  0r(f)  is  the  source 
laser  phase,  and  A^f)  is  the  envelope.  The  polarization  angle  between 
the  two  beams  is  <p(t).  likewise,  the  local  oscillator  laser  field  is 

E,  =  A,  cos[27rf(f  +  0,(0]  (2.2) 

where  f )  is  the  local  frequency,  0,(t)  is  the  local  phase,  and  A,  is  the 

(constant)  local  envelope.  The  photodetector  responds  to  the  combined 
fields  of  {r(0  +  /(f)}  with  intensity 

/c(f)  =  ~  A?(t) 4a,2+  Ar(f)A,  cos[2*(fr  -f,)  +  6r(t)  +  0,(n]cos^(O  (2.T) 
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The  output  of  the  photodetector  therefore  has  an  intermediate 
frequency  term  at  The  term  at  fr  +  f,  has  been  discarded  in  (2.3). 

The  intermediate  frequency  power  at  the  photodetector  will  be 

Pc(t)  =  2 JPrf  cos[2n(fr  -f,)t  +  0r(f)  -  0,(f)]cos<Mt)  (2.4) 

A  noncoherent,  direct  detection  system  responds  only  to  the  power 
of  the  incoming  signal,  with  an  intensity  of 

IM(t)  =  |  A,2(f)[l  +  cos(4  nfrt  +  20r(t))]  (2.5) 

The  double  frequency  term  will  normally  be  outside  the  photodetector’s 
response.  The  power  at  the  photodetector  will  be 

P«it)  =  Pr  (2.6) 

To  determine  why  coherent  systems  offer  superior  performance, 

the  quantum  performance  of  a  photodefector  must  be  considered.  A 
photodetector  requires  a  certain  number  Np  of  photons  to  fall  on  the 

detector  to  produce  a  single  electron.  For  an  ideal  detector,  Nf)  =  1.  The 
reciprocal  of  Np  is  called  the  quantum  efficiency  of  the  detector,  and  is 
denoted  by  tj.  The  current  produced  by  Mphotons/scc  is  i  =  Mr\q,  where 
q  is  the  electron  charge.  Each  photon  of  light  has  energy  hf  where  h  is 
Planck's  constant,  and  f  is  the  light  frequency.  Therefore,  the  incident 
power  P  =  Mfh  and  i  -  qqP/hf .  As  the  number  of  photons  falling  on  the 
photodetector  becomes  small,  the  discrete  nature  of  the  photon  to 
electron  conversion  events  leads  to  a  fluctuation  in  the  output  current 
that  becomes  more  severe  as  the  number  of  photons  falling  on  the 
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detector  decreases.  This  is  a  noise  process,  and  is  called  shot  noise.  For 
a  signal  bandwidth  b ,  the  shot  noise  process  has  power  (DeLange,  1968) 

(2.7) 


hf 


The  signal  power  is  i2R,  so  the  signal  to  shot  noise  ratio  is 

5  _  i2R  t ]F 
N  ~  Ns  ~  2 hfb 


(2.8) 


The  above  analysis  is  simplified,  in  that  noise  components  resulting 
from  dark  currents  and  thermal  noise  are  not  included.  A  more 
complete  analysis  may  be  found  in  (Keiser,  1991). 

Comparing  (2.4),  (2.6),  and  (2.8)  reveals  the  superiority  of  coherent 
systems.  Since  the  local  oscillator  power  carries  the  signal  modulation, 
the  photodetector  may  work  at  much  higher  signal  to  noise  ratios. 
Currently,  coherent  receivers  achieve  at  least  10  dB  better  sensitivity  for 
a  given  SNR  than  noncoherent  receivers. 


D.  DEMODULATORS 

Coherent  heterodyne  lightwave  systems  are  further  classified  as  to 
the  processing  at  EF.  In  coherent  demodulation,  the  phase  of  the  local 
oscillator  (within  the  detector,  not  the  local  oscillator  laser)  is  matched 
precisely  to  that  of  the  IF  signal  via  a  phase  locking  technique.  A 
coherent  2FSK  system  is  depicted  in  Figure  2.5.  In  this  figure,  the 
incoming  and  local  oscillators  have  equal  phase  0(f),  and  the  input  signal 
is  contaminated  by  noise  n(f).  The  coherent  system  is  the  optimum 
approach  when  the  phase  of  the  source  and  local  lasers  is  known  and 
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stable.  Unfortunately,  stable  phase  is  difficult  and  costly  to  achieve. 

The  alternative,  which  costs  3  dB  in  theoretical  signal  to  noise  ratio,  is  to 
use  noncoherent  demodulation  (Figure  2.6).  In  this  scheme,  the  phase  of 
the  IF  signal  is  irrelevant,  since  both  inphase  and  quadrature 
components  of  the  detection  process  are  squared  and  added. 

Fluctuation  in  the  source  laser  and  laser  local  oscillator  phase  (phase 
noise)  causes  fluctuation  in  the  IF  phase  and  still  degrades  the 
performance,  but  the  phase  noise  does  not  of  itself  complicate  the 
detection  process.  The  noncoherent  2FSK  system  depicted  in  Figure  2.6 
is  the  system  studied  in  this  thesis. 


Figure  2.5  Coherent  2FSK  Demodulation 
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Figure  2.6  Noncoherent  2FSK  demodulator 


III.  LASER  PHASE  NOISE 


For  a  variety  of  practical  reasons,  the  preferred  devices  for 
commercial  and  military  applications  of  fiber  optic  communications  are 
semiconductor  injection  lasers  (laser  diodes).  It  is  the  relative  spectral 
purity  of  the  laser  diode  that  permits  the  employment  of  coherent 
optical  modulation  schemes.  The  spectral  purity  is  not  absolute, 
however.  These  devices  are  generally  affected  by  phase  noise,  which 
broadens  the  spectrum  of  the  laser  and  can  severely  degrade 
performance. 

Earlier  it  was  stated  that  coherent  receivers  enjoy  a  10  dB  sensitivity 
advantage  over  direct  detection  systems.  In  practice,  this  advantage  is 
only  realizable  when  the  effects  of  phase  noise  are  controlled. 

Laser  phase  noise  is  caused  by  random  spontaneous  emissions 
within  the  lasing  medium.  Each  emission  event  causes  a  random  jump 
in  the  phase.  This  is  a  random  walk  process,  as  time  progresses,  the 
value  of  the  phase  will  wander  away  from  where  it  started.  The  mean 
squared  deviation  of  the  phase  grows  linearly  with  time.  It  has  been 
shown  (Salz,  1985)  that  the  common  Lorentzian  shape  of  the  laser 
spectrum  is  due  to  this  phase  noise  process. 

Since  the  mean  time  between  the  random  jumps  is  infinitesimal,  the 
phase  noise  process  becomes  in  the  limit  a  Wiener  process  characterized 
by  a  zero  mean  white  Gaussian  noise  /i(t)  with  some  two  sided  power 
spectral  density  N0 
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0(f)  =  2*r£  ju(T)dr  (3.1) 

with 

E{d2(t)}  =  £|47r2£^i(T1)dT,£jLi(T2)dT2| 

=  Mo  Jo  N°5(t'  '  T>  )dr\dr2  (3.2) 

=  4n2N0t 

To  determine  a  value  for  N0,  consider  the  laser  emission  process 

5(f)  =  cos[2nf0t  +  0(f)  +  0]  (3.3) 


where  <p  is  a  uniformly  distributed  random  variable. 
The  power  spectral  density  of  s(f)  is  (Salz,  1983) 


Gif) 


1 

1 

J 

1 

4  ir2N0 

1  + 

ff  +  f0) 

2  + 

1  +  1 

ff-fo  f 

V  *N0  J 

<  ttN0  j 

(3.4) 


This  is  the  Lorentzian  line  shape  shown  in  Figure  3.1.  Denoting  the  half 
power  linewidth  as  J3,  equation  (3.4)  may  be  solved  for 

N0  =  £  (3.5) 

The  random  phase  variable  0(f)  is  therefore  a  zero  mean  white  Gaussian 
process  with  PSD  from  (3.2)  of  2nf5t. 
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IV.  UNCODED  SYSTEM  PERFORMANCE 


The  coherent  system  with  noncoherent  demodulation  of  Figure  4.1 
will  be  considered  to  determine  the  deleterious  effects  of  phase  noise. 

A.  NONCOHERENT  2FSK  ERROR  PERFORMANCE  FORMULATION 

The  input  to  the  demodulator  is  the  IF  output  from  the  optical 
heterodyne  receiver 

r(t)  =  A  cos[2  jzf0t  +  0(f)]  +  n(t )  (4.1) 

It  is  assumed  that  a  data  0  is  transmitted.  The  signal  is 
contaminated  by  additive  white  Gaussian  noise  with  variance  N0/2. 

From  (3.2)  and  (3.5),  the  phase  noise  0(f)  is  white  Gaussian  with  variance 
2n(it .  It  is  assumed  that  the  frequency  spread  between  f0  and  is  large 
enough  to  assure  orthogonality. 

The  decision  variables  are 

Y°  =  exp[j0(t)]dt  + 

y,  =  K  +J"„f 

where  the  n  are  independent  identically  distributed  white  Gaussian 
random  variables  with  variance  N0T/ 4.  If  the  phase  noise  integral  term 

is  defined  as 

A  ,T  2 

Y2  =  —jo  exp[j'0(f)]df  (4.3) 


(4.2) 
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Figure  4.1  Noncoherent  Demodulator 


then  the  probability  of  error  conditioned  on  Y  yields  the  standard 
noncoherent  detection  problem.  The  density  for  Y0  conditioned  on  Y  is 
the  noncentral  x2  distribution 


kfW  =  ^exp| 


-(yo  +  i'2) 

2  a2 


Yjy~o' 


(4.4) 


where  I0  is  the  modified  Bessel  function  of  zero  order,  and 
a2  =  N0T/ 4. The  density  for  F,  is  the  exponential  density 


(4.5) 
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with  a2  =  N0T/ 4.  The  error  probability  conditioned  on  ^becomes 


Pr{v,  >  y„|y}  =  iexp 

1  J-Y r 

=  -exp  - 

2  "  N0rJ 


When  the  signal  is  not  corrupted  by  phase  noise,  the  bit  error  probability 


1  J  1  £„ 

2  2N0 

1  J-a2t~ 

-exp - 

2  14JV0J 


where  £h/NyJ  is  the  signal  energy  to  noise  ratio. 

The  expression  for  Y  may  be  rewritten 

.  <> 

Y2  =  j  £  exp[j0(t)}lt 

A*T2  rl  r  K 

=  —j-  j0  exp[j0(xT)]dx| 

=  ~~  exp[j'Vrv/(x)]cfx| 


(4.8) 


by  making  the  change  of  variable  t  =  x7\  The  random  variable  vdx)  has 
£{v/(x)j  =  x,  and  y  =  2kPT,  where  p  is  the  combined  3  dB  linewidths  of 

the  source  and  local  lasers  and  T  is  the  bit  rate.  If  the  random  variable 


X{y)  is  defined  as 


X(y)  =  I  exP[j 
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the  expression  for  bit  error  conditioned  on  X(y)  is 

-AzTX(y) 


-  F  J 

[fexp 

- 

i 

fl 

- 

—  F  J 

^X(Y)  | 

H 

4JV0 

JA 

2  Nn 


X(y) 


(4.10) 


From  (4.9),  it  is  clear  that  0  <  X(y)  <  1.  This  property  indicates  that 
A'(y)  may  be  considered  as  a  power  fraction  due  to  phase  noise. 

B.  PHASE  NOISE  VISUALIZATION 


Significant  insight  into  the  phase  noise  process  may  be  obtained 
from  careful  consideration  of  (4.9).  It  is  intuitive  that  when  the  phase 
angle  does  not  wander  during  an  integration  period  represented  by 
0  <  x  <  1,  the  integration  is  fully  coherent,  and  therefore  A(0)  =  1.  From 
previous  discussions,  the  phase  angle  y[yy/{x)  is  a  zero  mean  white 
Gaussian  process  with  a  variance  2tt/37x  =  yx  that  increases  throughout 

the  integration  period.  The  probability  density  for  the  process 
Z  =  ylyy/ix)  is 


&[*(*)] 


^2n yx 


exri- 


zz 

lyx 


<  Z  <  x 


(4.11) 


A  plot  of  the  density  versus  phase  angle  and  time  for  a  generic 
value  of  yis  shown  in  Figure  4.2. 

The  phase  angles  of  Figure  4.2  are  shown  in  a  continuous  form.  In 
actuality,  since  phase  angle  is  a  Modulo-2^  process,  the  density  function 
becomes 
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Figure  4.2  ft[z{x)] 


&[*<*)] 


^'Irzyx 


exrt- 


(z  +  2nrif 
2yx 


-  k  <  z  <  n 


(4.12) 


The  density  "rolls  up"  into  the  form  shown  in  Figure  4.3.  Now  the  impact 
of  phase  noise  becomes  very  evident,  in  that  as  the  integration  period 
progresses,  the  density  for  the  phase  begins  to  approai  h  a  uniform 
density,  which  destroys  the  coherence  of  the  integrati  u  Figures  4.4 
through  4.10  show  these  "rolled  up"  densities  for  a  rang-  <4  y.  As  y 

increases,  the  transition  to  a  uniform  density  is  more  pronounced  and 
rapid,  leading  to  less  time  for  coherent  integration  and  lower  signal 
power. 
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Figure  4.3  fz[z(x)]  with  phase  angle  taken  Modulo-2^ 


Figure  4.4  f/[z(xj]  for  y  =  1/4.  The  contours  of  the  2-D  plot  are  at  x 
intervals  of  (.2, .4, .(>,.8, 1.0) 


Figure  4.6  fz[z(x)]  for  y  =  1. 
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Figure  4.9  fz[z{x)\  for  y  =  8 


Figure  4.10  fz[z(x)\  for  y  =  16. 

The  above  intuitive  interpretation  provides  a  means  to  establish  an 
expected  value  of  the  random  variable  X(y).  Employing  (4.1 1),  this 
function,  call  it  X£(y),  may  be  defined  as 


X  E  (y )  =  £  £  expljz  |  fz  ( z  )dzdx 


A  plot  of  XE(y)  is  shown  in  Figure  4.1 1.  This  figure  demonstrates  the 
degradation  in  signal  power  due  to  increasing  phase  noise. 


Figure  4.11  Af(y) 


C.  PHASE  NOISE  PDF  COMPUTATION 

A  more  traditional  and  complete  analysis  of  the  bit  error  probability 
for  noncoherent  2FSK  requires  the  determination  of  the  actual 
probability  density  for  X(y).  This  is  a  formidable  problem,  and  to  this 

point  in  time,  only  rough  approximations  to  the  density  have  been 
obtained. 

The  first  approach  to  solving  for  the  PDF  of  A  was  put  forward  by 
(Foschini,  Greenstein,  and  Vannucci,  1989)  and  consisted  of  expanding 
the  RV  Ain  a  series  expansion,  retaining  the  linear  /  term 

A,  =  1  -  yz'Dz  (4.14) 
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where  z  is  an  infinite  dimension  vector  of  zero  mean,  unit  variance 
Gaussian  random  variables,  and  D  =  diagjd,,},  =  l/(i/r)2.  Examination 

of  XL  reveals  that  it  can  take  on  negative  values,  and  since  X  represents 
a  power  fraction  (with  permissible  values  0  <  X(y)  <  1) ,  this  characteristic 
is  obviously  unsatisfactory. 

What  is  required  then  is  an  approximation  to  X,  call  it  XA,  that  is 
easy  to  work  with  analytically,  has  a  range  of  values  0  <  X A(y)  <  1  for  all 
nonnegative  values  of  y ,  and  closely  matches  the  behavior  of  the  actual 

variable  X.  Another  way  of  stating  the  last  requirement  is  that  the 
moments  of  the  approximation  XA  match  those  of  X.  Such  an 

approximation  was  formulated  in  (Aziz  oglu  and  Humblet,  1991),  and  is 
given  by 

XA{y )  =  exp(->zTDz)  (4.15) 

where  z,D  are  as  given  above.  The  moments  for  the  approximation  are 
given  as 

for  all  real  t  >  ~nlj'ly . 

Given  the  moments,  the  PDF  of  XA  may  be  determined  through  a 
Gaussian  quadrature  technique.  The  general  technique  is  due  to 
(Hamming,  1973),  and  it  is  believed  that  the  application  of  this  technique 
to  this  class  of  problem  is  new.  Consider  the  Gaussian  quadrature 
problem 


fit 

sinh  yflyt 


(4.1G) 
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fK(x)f(x)dx  =  Ywkf(xk) 

Ja  tt 

K(x)  >  0 
a  <x<b 

If  K(x)  represents  an  arbitrary  probability  density  function  of  the 
random  variable  x,  the  moments  of  x  are  given  by 

m,  -  £{*'}  =  jhK(x)xtdx 

N 

=  L  wkK 

k=\ 

=  |  w 

\sinhyj2yt 

where  for  an  N  dimensional  Gaussian  quadrature,  t  =  0,  1,  ...  2 N  - 
define  the  following  NLh  order  polynomial  with  roots  xk 

tt(x)=  (x  —  x,  )(x  —  x2)-  *  •  (x  —  XN  ) 


(4.1 7) 


=  +cN_lxN-'  +  •••  +c0 


CN  =  1 


Next  multiply  each  t-  moment  equation  by  cf,  giving 

N 

=  2>*c  0 

N 

=  2>*X*C« 

N 

m2c2  =  £wfcxfo 


(4.18) 


2N-1.  Now 


(4.19) 


1 

JV 

mN  =  £wkx^ 


(4.20) 
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and  add  them  all  together 


Z™^  \  +  mN  =  ZZwiXc< 

(=0  )  f,0  k^l 

=  t^Z^r 


k=  1  f=0 


=  Zw^K) 


Now  multiply  the  equation  for  the  (t  +  /)-  moment  by  c, 

s 

™,CQ  =  ZW*X*Co 

k=l 

N 

=  Zw^'c, 

fc=l 

N 

ynl+2c2  =  ^wX2c 2 


(4.21) 


mN*l  \CN-\  -  Z^X*  CN  I 


(4.22) 


=  Z^ 


and  add  them  together 


Z^rwC,  +mN+/  =  ZZ^k^r 


N  N 


=  ZwXZKc, 

k  =  1  f^O 

N 

=  Zwkx>(xk) 


(4.23) 


for  all  /  =  0,  1,  ...  N  - 1.  Observation  of  the  left  side  of  the  above 
equation  indicates  that  it  forms  a  system  of  N  equations  in  N  unknowns 


when  the  moments  m,  are  known.  In  matrix  notation,  the  system 
becomes 


'  Mo 

m, 

-  "V. ‘ 

’  C0  ' 

~ms 

mx 

m2 

= 

mN 

.CSA. 

rm 2N-1. 

(4.24) 


With  the  coefficients  of  the  polynomial  n(x)  in  hand,  the  roots  xk  can  be 
numerically  computed.  This  operation  will  not  be  considered  here. 
Recall  the  equation  for  the  moments,  repeated  here 

m r  =  £{xf|  =  j*K(x)x'dx 

= 


(4.25) 


This  last  equation  can  be  put  into  matrix  form  as 


1 

1 

...  1  • 

‘w,‘ 

’  "to  ‘ 

*2 

-  4 

w2 

= 

m, 

VN- 1 
_Al 

VN-1 

A2 

VN-1 
•  •  •  AlV 

; wn . 

(4.2G) 


Examination  of  the  integral  and  summation  equations  for  mr 
indicates  the  correspondence  to  the  first  order  of  wk  and  K{x)dx ,  such 
that 

wk  =  K(xk)Axk 


K(xk)  = 


(4.27) 


Ax, 


In  essence,  this  technique  generates  a  rectangular  approximation 
to  the  integral  of  K(x).  A  better  (smoother)  method  of  extracting  K(x) 
from  the  computed  wk  and  xk  is  to  do  a  spline  interpolation  on  the  data 


set,  and  then  adjust  the  wk  such  that  the  integral  of  the  spline 
approximation  equals  1,  since  any  valid  probability  density  function 
must  integrate  to  1. 

Implementation  of  this  technique  is  hampered  by  the  ill 
conditioning  of  the  matrix  equations  and  root  finding  operation. 
Computation  in  standard  16  place  arithmetic  limits  the  N  dimension  of 
the  problem  to  about  8.  The  use  of  a  symbolic  algebra  system  with 
arbitrary  precision  allows  a  larger  N  dimension  to  be  used.  The  results 
of  Figure  4.12  were  generated  using  the  program  Mathematica  with 
N  =  16  and  128  digits  of  precision.  Comparison  with  the  results  of 
(Azizoglu  and  Humblet,  1991)  show  equivalence.  With  suitable 
computing  resources,  the  solution  for  a  given  N  could  be  obtained 
symbolically. 
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Figure  4.12  Probability  Density  Function  of  XA(y) 


D.  PERFORMANCE  WITH  PHASE  NOISE 


The  density  obtained  for  XA{y)  in  conjunction  with  Pb  from  (4.10) 
yields  the  error  performance  for  noncoherently  demodulated  2FSK 

WX,A(x)d*  <4-28> 

Performance  curves  for  several  values  of  y  are  shown  in  Figure  4.13. 
The  values  for  SNR  represent  the  signal  SNR,  since  it  is  assumed  that  one 
data  bit  is  represented  by  one  signal  bit. 

The  impact  of  increasing  data  rate  can  be  explored  through  Figure 
4.13.  Starting  at  a  point  on  a  curve  that  represents  a  particular  system, 
increasing  the  data  rate  shifts  the  operating  point  to  a  curve  at  lower  y , 

at  a  lower  SNR  that  represents  the  reduced  energy  per  bit  that  results 
from  higher  data  rates.  For  example,  beginning  at  y  =  4,  a  10  s  error  rate 

corresponds  to  a  SNR  of  32  dB.  Doubling  the  data  rate  results  in  a  drop 
in  SNR  of  3  dB  to  29  dB,  on  the  y  =  2  curve.  This  point  corresponds  to  an 
improved  bit  error  rate  of  5x10  7.  Also  consider  starting  at  y  =  1/2, 10  9 

bit  error  rate  at  18  dB.  Doubling  the  data  rate  results  in  a  worse  bit  error 
rate  of  10"  on  the  y  =  1/4  curve  at  15  dB.  In  general,  for  relatively  large 
y ,  data  rate  increases  lead  to  better  error  vs.  SNR  performance,  while  the 
converse  holds  at  low  y . 


I.og(P„) 


V.  CODED  PERFORMANCE 


A  convolutional  code  is  generated  by  passing  an  information 
sequence  through  a  linear  shift  register,  and  combining  the  shift  register 
cells  algebraically  into  several  output  channels.  The  diagram  of  Figure 
5.1  shows  a  general  implementation  of  a  convolutional  coder. 

A.  CONVOLUTIONAL  CODER 

Most  generally,  a  sequence  of  k  information  bits  enters  the  decoder 
at  stage  1,  and  is  subsequently  shifted  through  the  register  in  k  bit 
blocks.  The  total  number  of  blocks,  including  the  input  block,  is  L. 

There  are  a  total  of  Lk  individual  cells.  An  actual  implementation  would 
require  (L  - 1  )k  memory  cells,  since  the  first  stage  is  available  directly 
from  the  coder  input.  The  possible  contents  of  the  (L  - 1  )k  coder 
memory  cells  are  termed  the  states  of  the  coder.  There  are  2a'l)k 
possible  states  in  this  coder.  The  input  and  memory  cell  contents 
(hereafter  termed  the  coder  cells)  are  algebraically  combined  via  the  n 
code  generators,  which  are  usually  described  as  generator  polynomials. 
More  generally,  the  generators  may  be  described  as  a  sequence  of  Lk  l’s 
and  0’s.  Where  the  generator  is  1  indicates  that  the  corresponding  coder 
cell  is  connected  to  the  output  cell  represented  by  this  generator. 
Mathematically,  the  output  of  each  generator  is  the  inner  product 


Figure  5.1  A  convolutional  coder 

Modulo-2  of  the  generator  and  coder  cells.  The  coder  maps  the  k  input 
bits  onto  n  output  bits.  The  code  rate  of  the  code  is 

R  =  -  (5.1) 

n 


A  convolutional  coder  is  completely  specified  by  the  parameters  L, 
k ,  n,  and  the  generators,  which  are  commonly  given  in  octal.  The 
mapping  of  the  octal  representation  to  binary  is  shown  in  Table  5.1. 


TABLE 

5.1 

Octal 

Binary 

0 

000 

1 

001 

2 

010 

5 

Oil 

4 

100 

A 

101 

(> 

110 

7 

111 

The  value  of  L  is  termed  the  constraint  length  of  the  code.  The 
value  of  Lk  provides  a  measure  of  the  complexity  of  the  coder. 

From  observation  of  the  coder  in  Figure  5.1,  it  is  obvious  that  bits 
shifted  out  of  the  shift  registers  have  no  effect  on  the  coder  output.  The 

constraint  length  therefore  puts  a  limit  on  the  length  of  the  coder's 
impulse  response  (the  response  due  to  the  input  sequence  {1,0,0, 0,...} 

with  initial  state  all-zero),  that  being  the  product  Lk.  The  impulse 
response  h  =  {hvl\,...,hLk\  of  the  coder  is  analogous  to  that  of  any  other 

linear  time  invariant  system.  For  any  input  sequence  i  =  {z',,z2 . i7],  the 

output  of  the  coder  will  be  the  convolution  Modulo-2  of  z  with  h. 

B.  REPRESENTATION  OF  CONVOLUTIONAL  CODES 

There  are  three  equivalent  methods  of  describing  a  convolutional 
code:  the  tree  diagram,  the  trellis  diagram,  and  the  state  diagram.  Each 
will  be  important  to  understanding  the  coder  and  describing  its 
performance.  The  following  example  is  taken  from  (Proakis,  1989). 
Consider  the  L-  3,  k-  i,  n~  3  coder  shown  in  Figure  5.2.  The 

generator  sequences  are  given  as 

g,  =  {1,0,0} 
g^  =  {1,04} 
g,  =  {144} 


(5.2) 


Figure  5.2  Coder  for  I  =  3,  k  =  1,  n  =  3  code 

In  octal,  these  generators  would  be  given  as  (4,5,7).  The  theory  of 

generators  frequently  leads  to  the  expression  of  the  generator  sequences 

as  polynomials,  in  which  case 

9M)=  1 

g2(x)=  l  +  x2  (5.3) 

g3(x)=  1  +  x  +  x2 

First  consider  the  tree  diagram  for  this  code  shown  in  Figure  5.3. 
This  diagram  was  constructed  by  starting  at  the  all-zero  coder  state,  and 
describes  the  outputs  and  resulting  coder  states  for  a  sequence  of 
inputs.  The  upper  branch  out  of  each  coder  state  represents  the 
transition  due  to  an  input  bit  1,  the  lower  branch  to  a  bit  0.  Starting 
from  an  all  zero  loading,  if  a  bit  1  enters  the  coder,  the  state  transitions 
to  10.  The  output  of  the  coder  is  111.  If  another  1  enters  the  coder,  the 
path  through  the  tree  follows  the  lower  branch  again  to  state  11,  and  the 
resulting  output  is  110.  At  the  input  to  each  node  in  the  figure,  the 
upper  number  represents  the  coder  output,  and  the  lower  number 
represents  the  coder  state.  The  letter  A  represents  the  state  00,  B 


represents  01,  C  represents  10,  and  D  represents  11.  The  upper  branch 
out  of  a  node  is  always  the  zero  input,  and  the  lower  branch  is  the  1 
input.  This  process  continues  for  each  input  bit,  resulting  in  an  ever 
expanding  diagram. 

Examination  of  the  tree  diagram  shows  that  the  pattern  of  states 
and  outputs  begins  to  repeat  after  the  third  branch.  Since  there  are  a 
finite  number  of  unique  combinations  of  paths  (combinations)  of  input, 
output,  and  coder  state,  the  identical  points  on  the  tree  diagram  may  be 
merged  together.  Referring  to  Figure  5.3,  the  boxed  grouping  of  states 
may  be  merged  into  a  single  representation.  When  this  is  performed  for 
all  repeating  paths,  the  trellis  diagram  results.  The  trellis  diagram  is 
most  useful  for  understanding  the  paths  of  received  sequences  through 
a  decoder.  Figure  5.4  shows  the  trellis  diagram  for  the  1  =  3,  k  =  1,  n  =  3 
code.  In  the  trelhs  diagram,  the  upper  branch  out  of  a  node  always 
corresponds  to  a  0  input  bit,  and  the  lower  branch  to  a  1  input  bit.  The 
annotation  next  to  each  branch  is  the  output  sequence. 

The  final  representation  of  interest  is  the  state  diagram.  The  state 
diagram  is  a  signal  flow  graph  showing  the  possible  states  and  the 
transitions  between  them.  The  state  diagram  of  Figure  5.5  shows 
transitions  due  to  a  1  bit  as  dashed  lines,  and  those  due  to  a  0  input  as 
solid  lines. 


39 


<KU 


000 

C  000 

a  ooo 

00° 

non 

A\m 

/ill 

~\i  1 1 

/An. 

/Am 

\ 

\ 

oiiX 

oiiX 

OllY 

B#  \ 

V 

/\U)0/ 

//\0°V 

A\Oo\// 

\  ooj/ 

\  ooi// 

\  OOj// 

\  ooiAA 

c* 

Ilk 

Ilk/ 

Ilk/ 

Ilk/ 

ow  \ 

oik  \ 

Oik  N\ 

D* 

101 

A 101 

101 

1 

Steady  State 

Figure  5.4  Trellis  Diagram  for  L  =  3,  k  =  1,  n  =  3  code 


The  example  depicted  is  for  the  1  =  3,  k  =  1,  n  =  3  code.  For  any  rate 
k/n,  constraint  length  L  code,  there  will  be  2(i ~m  possible  states,  each 
with  2k  branches  entering  the  state,  and  2k  branches  leaving  the  state 
(after  steady  state  is  reached  in  the  trellis). 

The  state  diagram  representation  provides  a  mechanism  to  develop 
a  code  transfer  function  using  standard  signal  flow  graph  techniques. 

The  transfer  function  of  the  code  provides  a  mechanism  to  determine 
the  distance  properties  of  the  code,  and  to  determine  the  error 
conditions  that  go  uncorrected  by  the  code.  Since  the  codes  being 
considered  are  linear  codes  (the  shift  register  and  Modulo-2  adders  of 
the  coder  are  linear  processes),  the  transfer  function  may  be  constructed 
from  the  zero  state  to  the  zero  state  without  any  loss  of  generality.  Each 
branch  of  the  state  diagram  may  be  labeled  with  a  gain  function  JNW,DW°, 
where  w,  is  the  weight  of  the  input  sequence  causing  the  transition,  and 
wa  is  the  weight  of  the  resulting  output  sequence.  The  /  term  serves  to 
keep  track  of  the  number  of  state  transitions.  After  labeling  each 
branch,  a  set  of  simultaneous  equations  can  be  constructed  and  solved 
for  the  overall  transfer  function.  Keeping  in  mind  the  explosion  of  states 
and  inputs  for  large  values  of  L  and  k,  it  can  be  seen  that  the 
determination  of  the  transfer  function  can  be  difficult.  The  Appendix 
contains  a  more  detailed  discussion  of  the  procedure,  as  well  as  a 
computer  program  for  use  with  a  symbolic  algebra  system  to  eliminate 
the  tedium  of  these  calculations.  The  state  diagram  with  branches 
annotated  according  to  the  above  convention  is  shown  in  Figure  5.6. 
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Figure  5.6  State  diagram  in  transfer  function  form  for  L  =  3,  k  =  1,  n  =  3 
code 


The  system  of  gain  functions  for  the  example  code  is  as  follows: 

Xb  =  JDXc+JDXd 
Xc  =  JND'X^  +  JNDXb 
XD  =  JND2XC  +  JND2Xd 
XAn  =  JD2Xb 


The  self  loop  at  state  A  is  ignored  in  this  formulation  since  it  contributes 

nothing  to  the  distance  properties  of  the  code.  The  system  of  equations 
is  solved  for  X.^  /Aa  giving 


T(D,N,J ) 


1  -  JND2  -  J2ND 2 


(5.5) 


When  this  transfer  function  is  expanded  in  Taylor  series,  the  transfer 
function  becomes 


T(D,N,J)  =  pNDe  +  J4N2D*  +  JSN3D'°  +  (5.6) 

This  form  reveals  the  distance  properties  of  every  erroneous, 
uncorrectable  path  through  the  decoder.  The  first  term  in  (5.6)  reveals 
that  the  minimum  distance  of  the  example  code  (the  minimum  exponent 
of  D),  is  6.  A  full  explanation  of  the  individual  components  of  each  term 
follows.  The  explanation  is  based  on  the  assumption  that  the  all  zero 
code  sequence  is  the  transmitted  sequence. 

/Term:  The  exponent  of  J  is  the  length  of  the  information  sequence 

with  a  coded  representation  that  departs  the  all-zero  state  and 
subsequently  remerges  with  it.  Examination  of  the  trellis 
diagram  of  Figure  5.4  confirms  that  for  the  example  code,  this 
path  length  is  3. 

NTerm:  The  exponent  of  N is  the  number  of  l’s  in  the  true  information 
sequence  for  the  erroneous  path.  For  the  example,  the 
erroneous  path  corresponds  to  the  true  information  sequence 
100,  as  can  be  determined  from  the  trellis. 

D  Term:  The  exponent  of  D  is  the  distance  of  the  received  coded  bit 

sequence  from  the  all  zero  sequence.  For  the  example  this  is  6. 
The  coded  sequence  for  the  example  code  is  111  001  Oil. 
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C.  VITERBI  OPTIMUM  DECODING  ALGORITHM 


At  this  point,  the  tools  required  to  discuss  the  optimum  decoding  of 
convolutional  codes  are  in  hand,  and  discussion  will  shift  to  this  topic. 
The  optimum  decoding  technique  is  called  the  Viterbi  decoding 
algorithm.  The  path  of  quickest  comprehension  is  by  example,  so  the 
previous  L-  3,  k  =  1,  n-  3  code  will  continue  to  serve  this  function. 

Consider  the  trellis  diagram  from  Figure  5.4  and  two  of  the  possible 
paths  through  it 

r0  =  {0,0, 0,0, 0,0, 0,0,0} 

Kj  =  {1TT,0,0T,0JU}  3-/) 

The  path  r0  is  the  topmost  path,  and  path  r,  is  the  path  that 
remerges  with  the  all  zero  path  after  3  information  bits.  One  method  of 
decoding  is  to  periodically  take  a  sequence  of  coded  bits,  and  determine 
the  distance  between  the  sequence  and  every  possible  trellis  path 
corresponding  to  the  sequence  length.  The  output  of  the  decoder  would 
be  the  information  sequence  corresponding  to  the  code  sequence  closest 
in  distance  to  the  received  sequence.  The  measure  of  distance  can  be 
accomplished  in  several  ways.  If  a  hard  decision  is  taken  for  each  bit  as 
to  it  being  a  1  or  0,  the  distance  metric  is  the  Hamming  distance.  The 
Hamming  distance  is  defined  for  binary  sequences  as  the  number  of 
places  in  which  the  two  sequences  differ.  The  sequences  r0  and  r,  have  a 

Hamming  distance  of  6.  If  the  received  bits  are  quantized  into  a  number 
of  levels  without  taking  a  firm  decision  as  to  each  bit  being  a  1  or  0,  the 
Euclidean  distance  between  the  quantized  bits  becomes  the  metric.  This 
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approach  is  called  soft  decision  decoding,  and  generally  yields  superior 
performance  to  that  of  hard  decision  decoding,  at  the  cost  of  greater 
system  complexity.  This  thesis  will  hereafter  concern  itself  with  hard 
decision  decoding  for  the  reasons  mentioned  in  the  introduction. 

Continuing  with  the  example,  consider  all  paths  that  branch  from 
the  state  at  which  rQ  and  r,  merge.  Two  sets  of  paths  will  be  created,  one 
with  r0  at  the  beginning,  and  one  with  r,  at  the  beginning.  The  metrics 
of  the  extended  paths  will  be  added  equally  to  the  metrics  of  r0  and  r,, 
preserving  the  difference  between  the  metrics  of  r0  and  rv  This  implies 

that  the  total  number  of  possible  paths  to  keep  track  of  can  be  split  in 
two  at  the  point  where  r0  and  rx  merge  by  eliminating  the  path  with  the 

higher  metric  up  to  that  point.  Therefore,  at  every  state  where  two 
branches  merge,  one  entering  branch  can  be  discarded  based  on  its 
metric  statistics  up  to  that  point. 

Now  consider  a  very  long  trellis  with  a  corresponding  received 
sequence.  The  rightmost  bits  are  the  most  recent  inputs  to  the  decoder. 
At  every  stage,  when  two  candidate  paths  merge  at  a  state,  only  one  path 
survives.  Eventually,  a  stage  will  come  when  only  one  path  back  to  the 
beginning  of  the  trellis  survives.  This  path  represents  the  maximum 
likelihood  estimate  of  the  actual  transmitted  data  and  would  be  the 
output  of  the  decoder.  Since  the  current  coded  data  input  is  further  to 
the  right  of  this  stage,  a  time  delay  for  the  decoding  process  is  implicit. 
Figure  5.7  shows  several  steps  through  the  procedure  for  the 
L  =  3,  k  =  1,  n  =  3  example  code.  The  received  sequence  for  this  example 
was  {101  000  100  000  000  000  000}.  For  the  case  when  there  is  a  tie 


46 


101  OIK)  100  000 


101  000  100  000  000 


101  000  100  000  000  000  101  000  100  000  000  000  000 

Figure  5.7  Decoding  process  for  L  =  3,  k  =  1,  n  =  3  code 

between  path  metrics  merging  at  a  state,  the  tie  is  resolved  in  favor  of 
the  uppermost  incoming  branch.  Each  node  in  Figure  5.7  is  labeled  with 
the  path  metric  up  to  that  point. 

Examination  of  Figure  5.7  shows  how  two  candidate  paths  quickly 

appear.  Eventually,  the  lower  branch  is  eliminated,  and  the  output 
sequence  becomes  {0,0, 0,0,0}.  The  number  of  steps  required  to  achieve  a 

single  survivor  is  a  random  variable.  In  practice,  this  number  of  steps  is 
limited  to  a  practical  number,  and  the  path  with  the  lowest  metric  at  the 
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limit  is  the  output  path.  Refering  to  Figure  5.7,  an  appropriate  limit 
would  be  5  steps. 

The  process  described  above  appears  to  be  costly  in  that  a  large 
number  of  possible  paths  and  computed  metrics  for  each  path  must  be 
kept  in  memory.  Specifically,  for  a  rate  k/n  constraint  length  L  code,  the 

decoder  must  keep  track  of  2a-l)*  surviving  paths  with  corresponding 
path  metrics.  At  each  stage,  since  there  are  2*  entering  branches  at  each 
state,  a  like  number  of  branch  metrics  must  be  computed  and  added  to 
the  appropriate  path  metrics.  These  paths  are  then  compared,  and  the 
smallest  one  kept.  All  of  these  mathematical  operations  and  the  memory 
required  to  store  the  possible  paths  has  limited  the  application  of  this 
technique  to  codes  of  relatively  short  values  of  Lk,  and  lower  data  rates 
than  is  achievable  through  more  conventional  block  coding  techniques. 
The  use  of  a  soft  decoding  scheme  increases  the  memory  required  to 
store  the  metrics,  depending  on  the  number  of  quantization  levels,  and 
the  computation  of  branch  metrics  requires  squaring  and  square  root 
operations  instead  of  simple  addition  operations,  so  even  more 
restrictions  apply  to  soft  decoding.  Nevertheless,  as  memory  becomes 
cheaper,  and  logic  speeds  higher,  the  gap  between  capabilities  of  block 
and  convolutional  codes  should  decrease. 

D.  CODED  PERFORMANCE  ANALYSIS 

When  hard  binary  decisions  are  taken  on  the  incoming  coded  data 
before  decoding,  the  communications  channel  is  called  a  binary 
symmetric  channel.  In  the  Viterbi  decoder  discussed  above,  errors  in  the 
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received  coded  data  manifest  themselves  in  incorrect  path  metrics.  As  in 
the  determination  of  the  code  transfer  function,  it  may  be  assumed  that 
the  all  zero  codeword  is  the  actually  transmitted  codeword  without  loss 
of  generality.  Recall  the  information  that  the  transfer  function  provides 
about  incorrect  paths  that  remerge  with  the  all  zero  path  Each  one  of 
these  (infinite  in  number)  paths  has  a  distance  d  associated  with  it.  In 
this  case,  these  distances  are  each  paths'  distance  from  the  all  zero  codf* 
sequence.  It  is  crucial  to  realize  that  an  error  in  a  received  digit  where 
the  correct  and  incorrect  paths  are  the  same  has  equal  effects  on  the 
metrics  of  both,  and  can  therefore  be  ignored.  It  is  this  feature  that 
permits  us  to  consider  no  more  than  d  errors  even  though  more  than 
this  number  could  occur  on  a  path.  If  a  particular  distance  d  is  an  odd 

number,  the  all  zero  path  will  be  selected  as  the  surviving  path  at  the 
particular  node  if  less  than  (d  +  I)/2  errors  occur  in  the  coded  sequence 

up  to  that  point.  If  (d  + 1)/2  is  an  even  number  the  two  paths  will  have 

equal  metrics,  and  the  probability  of  making  the  proper  choice  of 
survivor  is  1/2  when  this  number  of  errors  occurs.  The  probability  then 

of  the  incorrect  path  surviving  is  a  binomial  probability  function 


PM )  = 


k  (J.l)  V»\  K  > 


d  odd 


j 

l 

k-Jfi,  l 


Z\  d  \  k  / 1  \ti-k  \  (  d  ^  j-2 


(3.8) 


p"(l-p)  dcve  n 


An  upper  bound  can  be  found  for  PJd)  as 


PM)  <  2jp(l-p) 


(5.9) 
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The  value  of  p  in  the  above  equations  is  the  channel  error  probability. 

The  above  error  probability  is  that  for  a  single  path.  In  general, 
there  are  ad  possible  paths  that  merge  at  the  particular  node  of  the 

trellis.  There  also  is  an  infinite  number  of  possible  merging  distances 
after  the  initial  dfree .  These  paths  are  not  independent  of  each  other,  but 

the  assumption  that  they  are  yields  a  pessimistic  estimate  for  the 
combined  error  probability 

Pe<±adP2(d)  (5.10) 


j=j 


free 


This  technique  is  termed  the  Union  Bound  Approach.  It  ignores  the 
intersections  of  the  different  paths. 

The  values  of  ad  may  be  obtained  from  the  transfer  function 
T(D,1,1).  Substituting  the  approximation  for  P2(a)  into  the  equation  for 
path  error  probability  yields 


P  <  T(D,U) 


(5.11) 


The  statistic  of  greatest  importance  is  the  information  bit  error  rate, 
which  may  be  obtained  as  follows.  When  an  incorrect  path  is  chosen,  the 
information  bits  in  that  path  will  differ  from  the  transmitted 
information.  Since  the  all  zero  information  path  was  the  one 
transmitted,  the  number  of  errors  will  be  the  number  of  l’s  in  the 
incorrectly  chosen  path.  The  information  bit  error  probability  for  a 
single  path  is  the  product  of  the  path  error  and  the  number  of  l's  in  the 

erroneous  information  sequence.  Consider  the  transfer  function  in  the 
following  form  with  the  ad  representing  the  number  of  patln  of  weight 

d  with  information  weight  w 


(5.12) 


T(D,N,  1)=  £ 

d~dfreg  I 

=  ^adDd  NflJ} 

d-dfree 


The  mean  number  of  l’s  in  the  information  sequence  corresponding  to  a 
path  is  given  by  f(d).  The  information  bit  error  probability  for  all  paths 


of  a  given  weight  may  be  over  bounded  by  the  mean 

P>  <  f(d)Pe 

<  j_ad((d)P2id) 

d=dfrte 

Recognizing  that 

=  j^adDdf(d)NfiJ)i 


(5.13) 


(5.14) 


and  employing  the  bound  on  P2(d)  of  (5.9),  the  bit  error  probability 


becomes 


P„< 


dT(D,N) 


dN 


,V=I.D=2vlp(l-p) 


(5.15) 


This  expression  for  Ph  is  the  expression  that  will  be  used  to  compute  the 
performance  of  the  laser  phase  noise  contaminated  bmary  FSK  system 
described  earlier. 


E.  CODES  STUDIED 

Four  codes  will  be  considered.  The  first  code  is  the  R  =  1/2,  L  =  1 
code  with  generators  (in  octal)  of  (5,7)  (Odenwalder,  1970).  The  transfer 
functions  for  this  code  are 
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DSNJ3 


dT(D,N) 


dN 


T(D,N,J)  = 


Ds 


1  -  DN/  2  -  DNJ 


(20 -ir 

=  D5+4D6  +  12D7+32D8  +  80D9  +  192D,o  +  ••• 


(5.16) 


(5.17) 


The  next  code  is  the  1?  =  1/2,  1  =  4  code  with  generators  (15,17) 

(Odenwalder,  1970).  The  transfer  functions  are 

D6J4N(D  +  JN-D2JN) 


T(D,N,J)  = 


1  -  DJN  -  DJ2N  -  D3J3N  +  D2J3N 2  -  ■ 
D4J3N2  -  D2J4N 2  +  D4J4N2 


(5.18) 


dT(D,N) 


dN 


n 


2D6  -  -  2P8  -f  P9  +  D 

(D!+2D-lf 

=  2D6  +7D7  +  18D8  +  49D9  +  130D10  +  333D11  +  ••• 


(5.19) 


The  third  code  is  the  R  =  1/4,  1  =  3  code  with  generators  (5, 7, 7, 7) 

(Larsen,  1973)  and  transfer  functions 

J3N(DU  +  Dl0JN  -  D'4JN) 


T{D,J,N)  = 
3T(D,N)\ 


dN 


1  -  D3JN  -  D3J2N  -  D2J3N2  +  D°J3N 2 
2 D10+Dn  -  D13  -  2D14  +  D17 


(D6  -  2D3  -  D2  + 1)2 


(5.20) 


(5.21) 


=  2D'°  +  D"  +  4D12  +9D13  +8D'4  +  25D 


»15 


The  final  code  is  the  R  =  1/4,  L  =  4  code  with  generators  (13,15,15,17) 
(Larsen,  1973)  and  transfer  functions 


T(D,J,N) 


D'3J4N(  1  +  DJN  -  D3JN  +  /2N2  -  D4J2N2) 

1  -  D3JN  -  D3J2N  -  DSJ3N  +  D6J3N2  -  D*J3N2  ~D*J4N2+  ■■■ 
D*J4N2  -  D7J4N3  +  D"J4N3  -  DSJSN3  +  D9JSN3 
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(5.22) 


4D13  +2D14  -  6D16  -5D17  +2D19  +6D20  +  - 


dT(D,N) 


dN 


D21  -  2D23  -  2D24  +  D 


24 


}27 


N=1  (Du  +  D9-D7-2D5-2D3  +  l) 

=  4  D13  +  2D14  +  10D16  +  3D17  +  16D18  +  34D19 


(5.23) 


These  codes  are  of  relatively  short  constraint  length;  the  current 
standard  in  commercial  production  for  microwave  radio  use  is  a 
constraint  length  of  7.  The  intended  application  for  these  codes  is  at  a 
high  data  rate,  so  it  seems  prudent  to  choose  codes  that  offer  lower 
complexity  of  implementation. 


F.  CODED  PERFORMANCE 

The  computation  of  information  bit  error  proceeds  as  outlined 
above,  realizing  that  the  value  of  p  comes  from  the  development  for 

uncoded  phase  noisy  binary  FSK.  When  coding  is  employed,  since  the 
signal  bit  duration  T  decreases  by  the  factor  R,  the  signal  to  noise  ratio 
also  decreases  by  R.  Also,  and  most  importantly,  the  value  for  y 

decreases  by  R.  The  error  performance  curves  based  on  uncoded 
/values  of  1,  4,  and  16  are  shown  in  Figures  5.8  through  5.10.  To 

facilitate  direct  comparisons  with  the  uncoded  system,  the  SNR  values 
are  with  respect  to  the  uncoded  system.  The  actual  signal  SNR  is 
reduced  by  the  code  rate. 

The  traditional  additive  white  Gaussian  noise  results  for  coding 

show  that  code  constraint  length  is  the  dominant  factor  in  improved 
performance  since  larger  L  produces  a  greater  dfree.  The  results  shown  in 
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Figures  5.8  through  5.10  show  that  code  rate  is  the  dominant 
performance  factor  in  the  presence  of  phase  noise.  This  code  rate  effect 
comes  about  from  the  replacement  of  a  single  information  bit,  with  its 
mostly  non  coherently  integrating  latter  portion,  by  several  code  bits, 
which  integrate  more  coherently. 

Figures  5.8  through  5.10  show  an  extremely  large  coding  gain. 
Performance  of  this  quality  enables  the  realization  of  practical 
heterodyne  systems,  operating  at  50  -  100  Mbps  per  wavelength  channel 
employing  the  new  high  speed  Viterbi  (hard  decision)  decoders.  With 
current  laser  diode  sources  exhibiting  10  -  100  MHz  linewidths,  a  200 
Mbps  signal  rate  results  in  1  <  y  <  4.  It  should  be  stated  that  the  bounds 

employed  in  the  coded  bit  error  are  somewhat  loose  bounds,  resulting  in 
a  pessimistic  coded  error  rate. 
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Figure  5.8  Coded  Performance  y  =  1 
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Figure  5.9  Coded  Performance  y  -  4 
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VI.  CONCLUSION 


The  comments  made  in  the  introduction  about  the  tradeoffs  that 
arise  from  the  use  of  coding  can  now  be  considered  in  more  detail.  It 
seems  clear  that  in  the  absence  of  any  constraints  on  signal  rate  and 
decoder  performance,  coding  offers  superior  performance.  Given 
constraints,  the  choice  to  use  coding  may  become  less  clear-cut. 

If  higher  data  rates  are  not  required,  and  the  source  laser  and 
decoders  support  the  higher  signal  rate  coding  requires,  then  coding 
offers  significant  performance  gains  as  compared  to  a  system  that 
increases  the  signal  rate  without  using  the  increased  data  rate  available. 
For  the  uncoded  system,  this  is  equivalent  to  having  unused  data 
capacity  on  the  link.  Table  6.1  summarizes  the  required  signal  SNR  for 
given  bit  error  rates  and  y  values  employing  uncoded,  coded,  and 
uncoded  signal-rate-increased  systems.  The  y  values  are  referenced  to 
the  uncoded  system. 

A.  HIGH  DATA  RATES  VERSUS  LOW  DATA  RATES 

Experimental  coherent  systems  have  been  run  at  modulation  rates 
well  over  1  Gbps.  The  fastest  decoders  are  currently  limited  to  200 
Mbps.  It  will  be  assumed  that  the  combined  laser  linewidth  of  source 
and  local  lasers  is  100  MHz  or  less.  A  data  rate  of  200  Mbps  corresponds 
to  a  y  of  about  4.  The  performance  of  an  uncoded  system  is  therefore 
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TABLE  6.1.  SIGNAL  SNR  REQUIRED 

y  =  1 

1(T5  1(T7  1(L9 

16  19  23 

11  12  14 

8  9  11 

14  16  18 

13  I  15  1  17 

quite  bad.  In  this  case,  a  signal  rate  increase  is  the  only  option,  unless 
one  considers  multiplexing  coded  systems.  The  total  power  required  in 
the  multiplexed  systems  to  achieve  a  given  bit  error  rate  is  less  than  that 
required  in  the  increased  signal  rate  system,  but  at  significantly  greater 
complexity.  In  general,  the  choices  to  be  made  for  high  data  rate 
systems  are  not  easy  ones. 

For  low  data  rate  systems,  the  choice  does  appear  to  be  easy.  So 
long  as  a  decoder  for  the  data  rate  desired  is  available,  coding  offers  a 
solution  to  the  phase  noise  problem. 

The  critical  components  in  all  this  are  the  source  and  local  lasers.  It 
is  a  very  high  priority  within  the  communication  industry  to  develop  the 
perfect  laser  for  coherent  communications.  Work  is  proceeding  slowly 
toward  this  goal.  Semiconductor  lasers  offering  linewidths  less  than  10 
MHz,  and  modulation  rates  higher  than  500  MHz  are  currently 
advertised. 


_ System _ 

Uncoded 

1  =  3,  R  =  l]2  Code 
L  =  3,  R  =  1/4  Code 
2x  Signal  Rate 
4x  Signal  Rate 
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Given  a  perfect  laser,  coding  will  still  find  application  to  combat 
other  sources  of  noise  in  communication  systems.  The  original 
applications  of  coding  were  to  combat  additive  white  Gaussian  noise  in 
communication  channels.  As  long  as  economical  decoders  operating  at 
high  enough  data  rates  are  available,  coding  will  find  application. 


B.  COMPARISON  WITH  DIRECT  DETECTION 


Since  it  has  been  established  that  a  coded  coherent  system  is  at 
least  theoretically  able  to  offer  excellent  performance,  it  is  worthwhile  to 
consider  how  well  such  a  system  compares  to  the  traditional  direct 
detection  system.  The  first  step  is  to  consider  the  error  rate  versus  SNR 
of  the  direct  detection  system  (Keiser,  1991). 


P„  =  ^  erfc| 


fl 

£h ) 

(6.1) 


This  performance  is  plotted  in  Figure  6.1. 

The  direct  detection  error  performance  is  roughly  equivalent  to 
uncoded  coherent  2FSK  with  y  <  1.  When  comparing  with  the  coded 
2FSK  systems,  it  is  seen  that  the  R  =  1/4  codes  perform  as  well  or  better 
for  all  y  <  16,  and  the  R  =  1/2  codes  perform  as  well  or  better  for  all  y  <  4. 
It  appears  as  though  the  10+  dB  receiver  sensitivity  advantage  of 
coherent  systems  may  now  be  realized  in  the  face  of  phase  noise. 

Coding  may  be  employed  on  the  direct  detection  system,  but  the 
tremendous  gains  of  the  coherent  system  will  not  be  realized,  since  the 
main  factor  in  that  effect  is  the  effective  reduction  in  y  brought  about  by 
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system  are  likely  to  be  equivalent  to  the  well-known  AWGN  channel 
results  (about  6  dB),  so  the  phase  noise  contaminated  coded  system  still 
offers  superior  performance. 


Figure  6.1  Direct  Detection  Pb 


C.  AN  APPLICATION  AND  FUTURE  RESEARCH 

This  thesis  has  shown  that  good  performance  is  attainable  at  low 
data  rates  with  a  coded  coherent  binary  FSK  system.  An  attempt  has 
been  made  to  point  out  at  the  pertinent  places  some  of  the 
implementation  issues  facing  the  employment  of  such  systems.  A 
hypothetical  application  of  such  a  system  will  now  be  discussed. 

There  are  military  applications  of  fiberoptic  links  that  have  been 
conceived  not  so  much  for  the  high  data  rates  available,  but  for  the  light 
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weight,  high  tensile  strength,  and  RFI  immunity  aspects  of  fiber. 

Consider  a  link  between  a  low  value  (or  expendable)  asset  and  a  receiver 
station.  Consider  a  medium  range  system  (50  Km)  that  requires  a 
relatively  low  data  rate  (less  than  50  Mbps).  There  may  also  be  a  link 
from  the  receiver  station  to  the  low  value  asset  (LVA),  but  it  is  not  of 
interest  (it  would  make  use  of  the  same  fiberoptic  cable).  What  is 
required  in  the  LVA  is  a  simple,  inexpensive,  low  power  consumption 
optical  transmission  system.  The  cost  and  complexity  of  the  receiving 
station  are  of  secondary  importance. 

Before  deciding  whether  a  coherent  or  direct  detection  system 
should  be  used,  consider  the  following.  The  range  requirement  and 
desire  for  the  lowest  power  system  mandate  the  use  of  single  mode  fiber. 
Chromatic  dispersion  effects  also  need  to  be  minimized,  dictating  the 
use  of  a  laser  diode  source  as  opposed  to  an  LED  (Keiser,  1991,  pp.  326- 
327).  As  a  result,  the  two  major  cost  and  simplicity  factors  favoring  a 
direct  detection  system  are  no  longer  part  of  the  equation. 

The  choice  between  coherent  and  direct  detection  systems  now 
becomes  whether  the  complexity  and  cost  of  the  coherent  system  justify 
the  better  performance.  Earlier  it  was  stated  that  the  frequencies  of  the 
source  and  local  oscillator  lasers  fluctuate  with  temperature  and  small 
drive  current  fluctuations.  If  the  controllers  for  these  parameters  can  be 
removed  from  the  LVA,  then  most  of  the  complexity  problem  within  the 
LVA  disappears.  To  remove  these  components  may  require  a  special, 
adaptive  local  oscillator  laser  and  demodulator  section  within  the 
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receiving  station,  and  the  effects  on  performance  of  such  a  scheme 
requires  further  investigation. 

Polarization  mismatch  between  the  source  and  local  lasers  is  also  a 
factor.  In  general,  any  twisting  or  bending  of  a  fiber  will  alter  the 
polarization  of  the  signal.  Receivers  that  are  insensitive  to  polarization 
also  require  further  investigation. 

Current  single  mode  fibers  offer  attenuation's  as  low  as  0.3  dB/Km. 
This  low  attenuation,  coupled  with  the  sensitivity  advantages  of  coherent 
systems,  offers  tremendous  increases  in  the  operational  range  of  systems 
that  employ  this  concept  of  operation.  Given  that  the  hurdles  of  phase 
noise,  frequency  stability,  and  polarization  sensitivity  can  be  overcome, 
entirely  new  classes  of  weapon  and  sensor  systems  may  become  reality. 


APPENDIX 


The  following  Mathematica  codes  (Wolfram,  1992)  implement  the 
formulations  for: 

The  Gaussian  quadrature  method  for  the  determination  of  the 
phase  noise  pdf 
Code  transfer  functions. 

A.  GAUSSIAN  QUADRATURE  METHOD  FOR  PDF  DETERMINATION 
FROM  MOMENTS. 

The  first  section  defines  the  moments  and  structure  of  the  matrices 
involved. 

u [g_, i_] : =Sqrt [  Sqrt[2  g  i ) /Sinh [ Sqrt [ 2  g  i ] ] ] 
u[Cj_,  0  ]  :  —  I 

MA[g_, n_] : =Table [u ( g, j  +  k]  ,  { j , 0, n-1 } ,  { k, 0, n-1 ) ] ; 

MB  f g_, n_]  : =Table [-u [g,  k]  ,  ( k, n, 2 n-1 } ] ; 

X [ n_ ] :=Table [x^i, {i, 0,n>] ; 

MD [ xk_List ] :=Table[xk~i, { i , 0 , Length [ xk] - 1 } ] ; 

ME [ g_, n_] : =Table [u [g, k] , { k, 0, n-1 ) ] ; 

The  next  section  computes  the  xk  and  wk 


prec=128 ; 
n=l  6/ 


.  .  .  for  1 28  digits  of  precision 
.  .  .  for  N  =  16 


g=  1 ;  ...  for  y  =  1 

A=N [MA [g, n] , prec] ; . . .  Remove  the  n  [•.prec]  statements  for  symbolic 

answers 

B=N [MB [ g, n] , prec] ; 

CC=LinearSolve [A, B] ; 

AppendTo [CC, 1 ] ; 
poly=CC  .  X ( n ] ; 

XK=x  /.  NSolve  [poly==0/ x,  prec]  ;  .  .  .  use  Solve  for  symbolic 

answers 

DD=N [MD [XK] , prec] ; 

EE=N [ME  [  g,  n] , prec] ; 

WK=LinearSolve [DD, EE] ; 

F=Transpose [ {XK, WK} ] ; 

The  data  generated  by  the  code  above  is  next  exponentially 
interpolated  (giving  straight  lines  on  a  log  plot),  and  an 
interpolatingFunction  object  is  obtained.  This  object  is  numerically 
integrated  to  find  its  norm,  and  finally  divided  by  its  norm,  so  that  it 
integrates  to  1. 

FI =F  /.  {a_,b_}->(a,Log[b]}; 

Fi= Interpol at ion [ FI [g] , Interpol a tionOrder->l ] ; 
norm=NIntegrate [Exp [ Fi [x] ]  ,  {x, 0, 1 } ] ; 


F[g,  x_]  :=Exp  [Fi  [x]  ]  /norm;  ...  The  final  PDF 
Plot  [F[g,  x]  ,  {x,  0,  1 }  ]  ;  ...  Plot  PDF 

B.  DETERMINATION  OF  CODE  TRANSFER  FUNCTIONS 
Algorithm 

1.  Separate  the  all  zero  state  into  two  states,  xt  and  xf,  such  that 
transitions  to  the  all  zero  state  are  attached  to  the  xf  state,  and  all 
transitions  from  the  all  zero  state  are  attached  to  x,. 

2.  For  every  possible  state  (there  will  be  2U  1)k  of  them) ,  determine  the 
coder  output  sequence  that  occurs  for  every  possible  input  sequence 
of  length  k  (there  will  be  2*  of  them).  Recall  that  the  coder  output 
sequence  is  the  dot  product  of  the  generator  sequence  with  the 
length  Lk  input  and  shift  register  sequence  (mod- 2).  Also  determine 
the  resulting  state  for  the  given  input  sequence. 

3.  Connect  the  original  state  to  the  resulting  state  with  an  arrow,  and 
annotate  the  connection  by  the  gain  functional  LNW,DW°  where  w,  is 
the  weight  of  the  input  sequence,  and  wo  is  the  weight  of  the 
corresponding  output  sequence. 

4.  Keep  track  of  these  transitions  and  path  gains,  and  form  a  system  of 
equations  that  relate  each  state  as  a  sum  of  other  states  with 
corresponding  gains.  For  given  L  and  k ,  there  will  be  2U  Uk  equations 

in  2a,,k  + 1  unknowns,  since  the  all  zero  state  was  split  into  the  states 
xf  and  xr  Each  equation  should  be  in  the  form 


Xj  =  +  gj2Xj2  +  •••  +  flf  .2kx  *  where  the  subscripted  states  and 

gains  refer  to  other  states.  In  other  words,  each  equation  has  2k 
components. 

5.  Solve  the  system  of  equations  by  eliminating  all  states  except  for  x, 
and  xf,  and  solve  for  the  ratio  xf/xf . 

It  is  obvious  that  the  task  of  finding  the  transfer  function  for  typical 
codes  is  impossible  without  the  aid  of  a  computer.  Even  with  a  symbolic 
math  system,  the  problem  is  daunting  for  large  values  of  Lk.  The  above 
algorithm  has  been  implemented  as  a  Mathematica  program,  and  used  to 
determine  the  code  transfer  functions  employed  in  this  research.  The 
code  is  listed  below. 

The  code  is  given  as  a  Mathematica  Function,  the  assumed  inputs 
are  the  values  of  L,  k,  and  the  generators  in  the  form  of  a  matrix  with  the 
rows  representing  the  generators  in  binary  form. 

BinaryConvCodeTransFunc [L_Integer, K_Integer, 

generators_?Matr ixQ] := 

Module [ { } , 

(*  L  is  constraint  length  *) 

(*  K  is  input  block  length  *) 
n=Length [generators ] ; 

(*  n  is  number  of  generators  *) 
numstates=2^ (K*L-K) ; 
ls=numstates+ 1 ; 

M=2/'K; 
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xsol=Table [ { 0 } , { Is  }  ] ; 
allstates=Table [ 

Take [ Join i Table [0, {K* (L-l ) } ] , 
IntegerDigits [i, 2]  ] , 

-K* (L-l) ] , 

{ i, 0, numstates-1 } ] ; 

allinputs=Table [ 

Take [Join [Table [0, {K} ] , 

IntegerDigits [ i, 2] ] , 

-K], 

{  i , 0 , M- 1 } ] ; 

For  [m=l,m<=M,m++, 
input=allinputs [ [m] ] ; 

For [ i=l, i<=numstates, i++, 
state=allstates [ [i] ] ; 
output=Table [0, { n } ] ; 

bigstate=Flatten [Prepend [state, input ] ] ; 
newstate=Take [bigstate, K* (L-l) ] ; 

For [ j=l, j<=n, j++, 
output [ [ j ] ] = 

Mod [Apply [Plus , bigstate* generators [ [ j ] ] ] , 2] ; 

1  ; 

wi=Apply [Plus, input] ; 
wo=Apply [ Plus, output ] ; 

{ { k } }=Position [allstates, newstate]  ; 
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if [k!=l, 

AppendTo [xsol [[ k} 3 , JJ  NNAwi  DD^wo  x[i]  ]]; 

If [k==l  &&  i ! =1 , 

AppendTo [xsol [[ Is ] ] # JJ  NNAwi  DDAwo  x[i]  ]]; 

] ; 

3  ; 

xsol=Apply [Plus, xsol, { 1 } ] ; 
solnset^Table [x [i 3  ==xsol [[i]],{i/2,ls}]; 
elims=Table [x [i ] , { i, 2, n urns fates } 3 ; 

Cancel  [  (x [ Is ]  /.  Solve [ solnset, x [ Is ], elims ]) /x [ 1 ] ] 
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